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Abstract 

Self-avoiding polymers in strictly two-dimensional (d = 2) melts are investigated by means of 
molecular dynamics simulation of a standard bead-spring model with chain lengths ranging up to 
N = 2048. The chains adopt compact configurations of typical size R(N) ~ N u with v = 1/d. 
The precise measurement of various distributions of internal chain distances allows a direct test 
of the contact exponents 0o = 3/8, Q\ = 1/2 and 02 = 3/4 predicted by Duplantier. Due to the 
segregation of the chains the ratio of end-to-end distance R e (N) and gyration radius R g (N) becomes 
R*(N)/R 2 g (N) « 5.3 < 6 for N » 100 and the chains are more spherical than Gaussian phantom 
chains. The second Legendre polynomial P2(s) of the bond vectors decays as P2(s) ~ l/s 1+i/ ® 2 
measuring thus the return probability of the chain after s steps. The irregular chain contours 
are shown to be characterized by a perimeter length L(N) ~ R(N) dp of fractal line dimension 
dp = d — O2 = 5/4. In agreement with the generalized Porod scattering of compact objects with 
fractal contour the Kratky representation of the intramolecular structure factor F(q) reveals a 
strong non-monotonous behavior with q d F(q) ~ l/(qR(N))® 2 in the intermediate regime of the 
wave vector q. This may allow to confirm the predicted contour fractality in a real experiment. 
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I. INTRODUCTION 



Dense self-avoiding polymers in two dimensions (2D) have been considered theoretically 
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17 



by means of computer simulation 5 
2 j |. It is now generally accepted {2-6 



.0 



Jl6| and more recently even in real experiments 



11 



16 



17 



of Refs. HI and 



201 ] (with the notable exception 



18) ) that these chains adopt compact and segregated conformations at high 



densities, i.e. as first suggested by de Gennes [l( the typical chain size R(N) scales as 



N « pR d (N) 



(1) 



with N being the chain length, p the monomer number density and d = 2 the spatial dimen- 
sion. We assume here that monomer overlap and chain intersections are strictly forbidden 



[3] 



Compactness and segregation are expected to apply not only on the scale of the total 



n 



chain of N monomers but also to subchains comprising s monomers |3|, |20| , at least as long 
as s is not too small. The typical size R(s) of subchains should thus scale as 



R(s) « {s/pY for 1< s < N 



(2) 



with a Flory exponent v — l/d — 1/2 set by the spatial dimension. Interestingly, the direct 



visualization of chain conformations is possible for DNA molecules 17, 23, 241, nanorope 

n n 

polymer chains [18J or brushlike polymers [20[ adsorbed on strongly attractive surfaces or 



confined in thin films by means of fluorescence microscopy 17|, |23| or atomic force microscopy 



18 



20J. The experimental verification of the various conformational properties discussed in 



this paper, such the one described by Eq. (j21), is thus conceivable for these systems 



20] 



Compactness does obviously not imply Gaussian chain statistics [l| (as incorrectly stated, 



e.g., m 



181 ] ) nor does segregation of chains and subchains impose disk-like shapes minimizing 



the average perimeter length L(s) of chains and subchains (3). The contour boundaries are 
in fact found to be highly irregular as is revealed by the snapshots presented in Fig. [T] and 
Fig. |2J Elaborating a short communication made recently 16) we present here theoretical 
arguments and molecular dynamics simulations demonstrating that the contours are in fact 
fractal, scaling as 

L(s) ~ R(s) d * ~ s l - u@ \ (3) 

with dp = d — 62 = 5/4 > 1 being the fractal line dimension [25]. Our work is based on the 
pioneering work by Duplantier who predicted a contact exponent 2 = 3/4 [2] and a more 



recent paper by Semenov and one of us (A.J.) [3]. In contrast to many other possibilities to 
characterize numerically the compact chain conformations the perimeter length is of interest 
since it can be related to the intrachain structure factor F(q) with q being the wave vector. 
It is thus accessible experimentally, at least in principle, by means of small-angle neutron 



scattering experiments to all polymers which can be appropriately labeled 26J|. Specifically, 



it will be demonstrated that due to the generalized Porod scattering of compact objects 



26 



the structure factor of dense 2D polymers should scale in the intermediate wave vector 
regime as 

F(q) « N/(qR(N)) d+@2 ~ N^q- 11 ^ (4) 
x>r sufficiently long chains and not as F(q) ~ N°q~ 2 as numerous authors have assumed 

, ti- 



ll 



The present paper is organized as follows. In Sec. [TT] we recall the computational model 
used for this study. Our mainly numerical results are presented in Sec. IIHI We confirm 
first the compactness of the chain conformations by considering the typical size of chains 
and subchains (Sec. IIII A|) . Corrections to this leading power-law behavior due to chain-end 
effects caused by the confinement of the chains will be analyzed in Sec. IIII CI That v = 1/2 
does not imply Gaussian chain statistics will be emphasized by the scaling analysis of various 
intrachain properties such as the histograms of inner chain distances Gi(r) in Sec. IIII the 
bond-bond correlation functions Pi(s) and ^(s) i n Sec. IIII Dl or the single chain structure 
factor F(q) in Sec. IIII Gl Two (related) scaling arguments will be given in Sec. IIII Fl and at 
the end of Sec. IIII Gl demonstrating the scaling of the perimeter length L(s), Eq. (jHJ). The 
analytic calculation of the structure factor F(q) for 2D melts is relegated to the Appendix. 
A discussion of possible consequences of the observed static properties for the dynamics of 
dense 2D solutions and melts concludes the paper in Sec. HVl 

II. COMPUTATIONAL DETAILS 

Our numerical results are obtained by standard molecular dynamics simulations of 
monodisperse linear chains at high densities. The coarse-grained polymer model Hamilto- 
nian is essentially identical to the standard Kremer-Grest (KG) bead-spring model 0, 3(3] 
which has been used in numerous simulation studies of diverse problems in polymer physics 



29M32| . The non-bonded excluded volume interactions are represented by a purely repulsive 



3 



Lennard- Jones potential 33] 



U nh (r) = 4e [{o/r) 12 - [a/rf] + e for r/a < 2 1 / 6 (5) 

and U n b(r) = elsewhere. The Lennard- Jones potential does not act between adjacent 
monomers of a chain which are topologically connected by a simple harmonic spring potential 

U h {r) = k h (r - l h ) 2 (6) 

with a spring constant = 338e and a bond reference length = 0.967a. Both constants 
have been calibrated to the "finite extendible nonlinear elastic" (FENE) springs of the 
original KG model. Lennard- Jones units are used throughout this paper (e = a = 1). 
The classical equations of motion of the multichain system are solved via the Velocity- 
Verlet algorithm at constant temperature using a Langevin thermostat with friction constant 



7 = 0.5 



33] 



We focus in this presentation on melts of density p = 7/8 = 0.875 at temperature T = 1. 
Due to the excluded volume potential monomer overlap is strongly reduced as may be seen 
from the pair correlation function g(r) shown in Fig. [3l The bonding potential, Eq. ([6]), 
prevents the long range correlations which would otherwise occur at such a high density 



for a 2D system of monodisperse Lennard- Jones beads [331 . |34|. As shown in the inset of 
Fig. [3] we have in fact a dense liquid and the oscillations of g(r) — 1 decay rapidly with 
an exponential cut-off. The parameters of the model Hamiltonian and the chosen density 
and temperature makes chain intersections impossible, as can be seen from the snapshot of 
"chain 1" presented in Fig. [TJ We simulate thus "self-avoiding walks" in the sense of the 
first model class discussed in Ref. [3J]. 

Monodisperse systems with chain lengths N ranging between N = 32 up to iV = 2048 
have been sampled using periodic square boxes containing either 98304 or 196608 monomers. 
The larger box of linear length 474.02 corresponds to 96 chains of length N = 2048. Some 
conformational properties discussed below are summarized in the Table. Except the systems 
with iV = 2048, all chains have diffused over at least 10 times their radius of gyration R g {N) 
providing thus sufficiently good statistics. Note that our largest chain is about an order of 
magnitude larger than the largest chains used in previous computational studies of dense 
2D melts: N = 59 by Baumgartner in 1982 iV = 100 by Carmesin and Kremer in their 
seminal work in 1990 6], N = 100 by Nelson et al. in 1997 3, N = 32 by Polanowski 



and Pakula in 2002 Q, N = 60 by Balabaev et al. in 2002 |l0j, A" = 256 by Yethiraj in 
2003 Q and iV = 256 by Cavallo et al. in 2005 Q, Q- The presented data was obtained 



on IBM power 6 with the LAMMPS Version 21May2008 [35J. It is part of a broader study 
where we have systematically varied density, system size and friction coefficient to confirm 
the robustness of theory and simulation with respect to these parameters. Since the chain 
length is computationally the limiting factor fixing the number of "blobs" jlj at a given 
density we present data at the largest density, i.e. the largest number of blobs, where we 
have been able to equilibrate chains of AT = 2048. 

III. NUMERICAL RESULTS 

A. Chain and subchain size: Compactness 

Figure H] confirms that 2D chains are indeed compact as stated in Eq. (TTJ and Eq. (|2]) and 
as shown in Fig. [Tj for chains and in Fig. [2] for subchains of arc-length s. The typical size 
R(s) of subchains is characterized by either the root-mean-square end-to-end distance R e (s) 
or the radius of gyration R g (s) jsg], As indicated in the sketch we consider a subchain 
between two monomers n and m = n + s — 1 and average over all pairs (n, m) possible in a 



chain of length A" following [37H39I]. Averaging only over subchains at the curvilinear chain 
center (n,m « N/2) slightly reduces chain end effects, however the difference is negligible 
for the larger chains, A" > 1024, we focus on. The limit s = N corresponds obviously 
to the standard end-to-end vector R e {N) and gyration radius R g {N) of the total chain. 
Open symbols refer to subchains of length s < N with A" = 1024 (squares) and A" = 2048 
spheres), stars to total chain properties (s = N). In agreement with various numerical 
14] and experimental studies [l3, H] the presented data confirms that the chains 




5|,16J 



11 



are compact, i.e. v = 1/d (thin lines), and this on all length scales with 1 s < N. 

The segregation of the chains may also be shown by computing the average number n c h a in 
of chains in contact with a reference chain, i.e. having at least one monomer closer than 
a distance a 2 to a monomer of the reference chain (Sec. IIIIFI) . Depending weakly on 
a, we find n c h a in ~ 6, as one may expect for 2D colloids and in agreement with Fig. [TJ At 
variance to any open non-segregated polymer-like structure n c hain does thus not increase 
with chain length N. An alternative way to confirm this statement is to count the centers 



5 



of mass around the reference chain's center of mass by integrating the center-of-mass pair- 
correlation function <7cm(r), a structureless function without oscillations (not shown). One 
verifies that 6 chains are found in a shell of about 2R g (N) and this irrespective of N. 



B. Segment size distributions and contact exponents 

Being characterized by the same Flory exponent v as their three dimensional (3D) coun- 
terparts does by no means imply that 2D melts are Gaussian 0, Q]. This can be directly 
seen, e.g., from the different probability distributions of the intrachain vectors r = r m — r n 
presented in Fig. To simplify the plot we have focused on the two longest chains N = 1024 
and N = 2048 we have simulated. As illustrated in panel (a), 

• Go(r, s = N) characterizes the distribution of the total chain end-to-end vector (n = 1, 
m = N), 

• G*i(r, s = N/2) the distance between a chain end and a monomer in the middle of the 
chain (n = 1, m = N/2), 

• G2{r,s = N/2) the distribution of an inner segment vector between the monomers 
n = N/4 and m = 3JV/4, 

• while the "segmental size distribution" G e (r,s) [37] averages over all pairs of monomers 
(n, m = n + s — 1) for s < N. 

The second index s indicated in Gi(r,s) characterizes the length of the subchain between 
the two monomers n and m. As shown in panel (b), all data for different iV and s collapse 
on three distinct master curves if the axes are made dimensionless using the second moment 
R? of the respective distribution. The only relevant length scale is thus the typical size 
of the subchain itself. The distributions are all non-monotonous and are thus qualitatively 
different from the Gaussian (thin line) expected for uncorrelated ideal chains. Confirming 
Duplantier's predictions [2!] we find 

RtG l (r,s) = x *f l (x) (7) 

with x = r/Ri being the scaling variable and the contact exponents 60 = 3/8, 61 = 1/2 
and 02 = 3/4 (dashed lines) describing the small-x limit where the universal functions 



fi(x) become constant. Especially the largest of these exponents, G2, is clearly visible. The 
contact probability for two monomers of a chain in a 2D melt is thus strongly suppressed 
compared to ideal chain statistics (6 = ©i = 6 2 = 0). 

The rescaled distributions show exponential cut-offs for large distances. The Redner-des 
Cloizeaux formula 40] is a useful interpolating formula which supposes that 

fi(x) = Ciexp(-kiX 2 ). (8) 

The constants hi = 1 + Qj/2 and q = /irYiJzi) with Y{z) being the Gamma function 41] 



are imposed by the normalization and the second moment of the distributions |42j. This 
formula is by no means rigorous but yields reasonable parameter free fits as it is shown by 
the solid line for f2(x). 

Obviously, G e (r,s) ~ G (r,N) for very large subchains s N (not shown). As can be 
seen, the rescaled distributions G e (r, s) and G2(r, N/2) become identical if 1 <C s <C N. It is 
for this reason that the exponent B2 is central for asymptotically long chains as will become 
obvious below in Sec. IIIIFI and Sec. IIII Gl The two exponents Go and Gi are only relevant 
if the measured property specifically highlightes chain end effects as in the example given in 
the next subsection. 



C. Chain and subchain size: Corrections to asymptotic scaling 

The log-log representation chosen in Fig. H] masks deliberately small corrections to the 
leading power law due to chain end effects which exist in 2D as they do in 3D melts {3?]]. 
These are revealed in Fig. [6] presenting R e (s)/s 1 ^ 2 and R g (s)/s 1 ^ 2 vs. s and R e (N)/N l l 2 
and R g (N)/N 1 / 2 vs. N using log-linear coordinates and the same symbols as in Fig. 0} The 
reduced radius of gyration (bottom data) becomes in fact rapidly constant and chain length 
independent. As emphasized by the bottom horizontal line we find 

b g « R g (s)/s 1/2 « R g (N)/N 1/2 for s, N > 100 (9) 

with b g = 0.65 for p = 7/8. Interestingly, we observe non-monotonous behavior for 
R e (s)/s 1 ' 2 with a decay for s > N/2. Due to this decay R e (N)/N 1 ^ 2 (stars) is system- 



atically below the corresponding internal chain distance R ( 



s) / s 1 / 2 . If fitted from the chain 



end-to-end distances one obtains an effective segment size 36| 



b e * w R e {N)/N 1/2 w 1.5 for 100 < N, (10) 



111. If 



as shown by the dashed line. (The index * indicates that we refer to the total chain.) This 
value corresponds to the ratio Rl(N)/R 2 (N) ~ 5.3 < 6 given in the Table. It confirms 
similar observations made in previous simulations using much shorter chains 5|, |6j, [8 
on the other hand the effective segment size b e is obtained from the internal distances this 
yields 

b c « R e (s)/s l/2 » 1.6 for 100 < s < N, (11) 
as indicated by the top solid line. This value is consistent with a ratio R 2 (s)/R 2 (s) = 



(b e /b g ) 2 ~ 6 as one expects in any dimension due to 36| 

R 2 g (s) = I f ds\s - s')Rl(s') = h 2 e s (12) 
s Jo 

where we have assumed R 2 (s') = b 2 s' for all s' up to s' = s. Since this assumption breaks 

down for s — > N (as seen in Fig. [6]) Eq. f|T2|) is not in conflict with Eq. ([TO]) . Note that 

the integral over s' in Eq. (|T2|) is dominated by subchains with s' ~ s/2 and that thus 

large subchains of order s' ~ N are less relevant for the gyration radius of the total chain 

R g {s = N). Hence, the non-monotonous behavior observed for R e (s) should be barely 

detectible for the radius of gyration in agreement with Eq. ([9]). 

That a naive fit of b c from the total chain end-to-end distance R e {N) leads to a systematic 

underestimation of the effective segment size of asymptotically long chains is a well-known 



fact for 3D melts 37 1. However, both estimations of b e merge for 3D melts if sufficiently long 
chains are computed (as may be seen from Eq. (16) and Fig. 4 of Ref. [33]). Apparently, 
this is not the case in 2D since if R e {s) or R 2 (s)/s are plotted as a function of x = s/N a 
nice scaling collapse of the data is obtained for large x and for iV > 256, i.e. 

RI(S) » b 2 e S 2 »f(x) With f{x) = \ 1 lf X « 1 (13) 

(Since a very similar scaling plot is presented in the inset of Fig. [7] this figure is not given.) 
Hence, chain end effects may not scale away with N — > oo as they do in 3D. A simple 
qualitative explanation for Eq. ( ITBl is in fact readily given by considering an ideal chain of 
bond length b e squeezed into a more or less spherical container of size R c ~ iV 1 / 2 . For s <^ N 
it is unlikely that the subchain interacts with the container walls and R e (s) ~ b e s 1 ^ 2 . With 
s — >■ iV the chain will feel increasingly the confinement reflecting it back from the walls into 
the center of the container reducing thus the effective segment length b eir associated with 

8 



the chain end-to-end distance. Since the scaling function f(x) must be universal, it should 
be possible to express the ratio {b c +/b e ) 2 — and thus the ratio Rl(N)/R 2 (N) = (b eit /b g ) 2 
- in terms of the dimension d and universal compact exponents Go, 0i and 02- The two 
statistical segment sizes b e and b eic should thus be related. At present we are still lacking a 
solid theoretical proof for the latter conjection. 



D. Intrachain orientational correlations 



Let e$ denote the normalized tangent vector connecting the monomers i and i + 1 of 
a chain. The bond-bond correlation function P\{s) = (e n • e m ) has been shown to be of 
particular interest for characterizing the deviations from Gaussianity in 3D polymer melts 



371 . |38| . (As above we average over all pairs of monomers (n, m) with m = n + s — 1.) The 



reason for this is that 



38] 

d 2 Rl(s 



#^ (14) 

and that thus small deviations from the asymptotic exponent 2v — 1 are emphasized j^ . 
P\{s) is presented in Fig. [7]for different chain lengths. Apparently, the bond pairs are aligned 
only for small arc-lengths with s < 10. For larger s the bonds are anti- correlated with two 
characteristic peaks. The first anti-correlation peak visible in panel (a) is due to the local 
backfolding of the chain contour which can be directly seen from chain 1 drawn in Fig. [TJ 
Note that the chain length dependence disappears for N > 256. The second anti-correlation 
peak is shown in panel (b) where —P\(s)N is plotted as function of x = s/N. Using Eq. ( |T4l) 
this corresponds exactly to the scaling expected from Eq. ( fT3l) with an associated universal 
function scaling as —d 2 (xf(x)). In agreement with the qualitative explanation mentioned 
at the end of Sec. IIII Al the peak at s ~ N/2 can be attributed to the confinement of the 
chain which causes long segments to be reflected back, i.e. Pi{s) must be anti-correlated 
for large x. (Data for iV = 2048 is not presented here due to its insufficient statistics.) We 
stress finally that altogether this is a rather small effect and essentially Pi(s) ~ d 2 R 2 (s) w 
for s ^> 10 2 and iV ^> 10 3 as already obvious from panel (a). Hence, the first Legendre 
polynomial confirms that to leading order 2u w 1 for sufficiently large chains and segments. 

Conceptually more important for the present study is the fact that the second Legendre 
polynomial ^(s) — (( e n ' e m) 2 ) — 1/2 given in the main panel of Fig. [8] reveals a clear power 
law behavior over two orders of magnitude in s (dashed line). This power law is due to 
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(a) the return probability p r (s) after s steps and (b) the "nematic alignment" of two near- 
by bonds. The alignment of bonds is investigated in the inset of Fig. [8] where the second 
Legendre polynomial is plotted as a function of the distance r between the mid-points of 
both bonds. Averages are taken over all intrachain bond pairs with [r, r + Sr] using a bin 
of width 5r = 0.01. Since i^tj) becomes rapidly chain length independent we only indicate 
data for N = 2048. The vertical axis is rescaled with the phase volume 1-nr. As can be seen, 
the orientational correlations oscillate with r and this with a rapidly decaying amplitude. 
These oscillations are related to the oscillations of the pair correlation function g(r) shown 
in Fig. Eland reflect the local packing and wrapping of chains composed of discrete spherical 
beads. Due to both the oscillations and the decay only bond pairs at r pa 1 matter if we 
compute P2{s), i.e. if we sum over all distances r at a fixed curvilinear distance s. Following 
Eq. ([7j) one thus expects 

Pi{s) ^Pr(s) = limG e (r,s) ~ l/s 1+u&2 = l/s 11/s (15) 

r— >0 

for 1 <C s <C N. The agreement of the data with this power law is excellent and provides, 
hence, an independent confirmation of the contact exponent 02 = 3/4. 

E. Chain and subchain shape 

As obvious from Fig. [I] and Fig. [2] the conformations of chains and subchains are neither 
perfectly spherical nor extremely elongated. Having discussed above the chain size we ad- 
dress now the chain shape as characterized by the average aspherity of the gyration tensor. 
The gyration tensor M of a subchain between the monomers n and m = n + s is given by 

^ m 

M a/3 = ~ y^X r i,a ~ R cm,a)(ri,(3 ~ #cm,/s) (16) 
i=n 

with R cmt a being the a-component of the subchains's center of mass. We remind that the 
radius of gyration Rg(s) discussed in Sec. IIII Al is given by the trace tr(M) = M xx + M yy = 
Ai + A2 averaged over all subchains and chains with eigenvalues Ai and A2 obtained from 

Ai, 2 = X - (tr(M) ± y/tr^) 2 - 4det(M)) . (17) 

The ratio of the mean eigenvalues (Ai) / (A2) for s = N is given in the Table. Decreasing 
slightly with this "aspect ratio" approaches 

(Ai) : (A 2 ) « 4.5 : 1 (18) 
10 



for our longest chains which corresponds to a reduced principal eigenvalue (Ai) 0.8. It 



should be noted that Gaussian chains and dilute good so 
by an aspect ratio (Ai) / (A 2 ) ~ 5.2 and 6.7, respectively 
elongated. The asphericity of t 
by computing the moments 24 



vent chains in 2D are characterized 



431 ] . Our chains are thus clearly less 



re inertia tensor of 2D objects may be further characterized 
43-45| 



(Ai - A 2 > 



, A 2 (s) 



((Ai - A 2 ) 2 > 



(19) 



(Ai + A 2 )' zw ((Ai + A 2 ) 2 ) 

which are plotted in Fig. for subchains (s < N) and total chains (s = N) using the 
same symbols as in previous plots. Ai = 2 (Ai ) /R 2 Z — 1 describes the mean ellipticity and 



44 



45| . Obviously, Ai = A 2 = 1 for rods and 



A 2 the normalized variance of Ai and A 2 
Ai = A 2 = for spheres. Note that taking the first and the second moments of the ellipticity 
(Ai — A 2 )/(Ai + A 2 ) of each subchain yields qualitatively similar results (not shown). As one 
expects, both moments do not depend on whether a chain or a subchain is considered. In 
agreement with Yethiraj Uj they decrease weakly with N and s. Unfortunately, it is difficult 
to determine precisely the plateau values one expects for asymptotically large chains and 
subchains, and the horizontal lines with 



Ai = 0.63 and A 2 = 0.51 



(20) 



are merely guides to the eye. Note that Yethiraj ll|] indicates A 2 = 0.52 for N = 256. 
Considering that the latter value has been obtained at a slightly smaller monomer volume 
fraction both A 2 -values are compatible. We remind that in two dimensions Ai ft 0.68 and 



and Ai ft 0.74 and A 2 ft 0.62-0.64 



43 



45] . These values are definitely 



A 2 = 2(d+2)/(5d+4) ft 0.57 for Gaussian chains 43j,|44 
for dilute good solvent chains according to Refs. Ill, [24 
larger than our respective estimates, Eq. (|20|) . and segregated chains in 2D melts are thus 
clearly more axisymmetric. 

The above analysis has been motivated by recent experimental work on the conforma- 
tional properties of dilute DNA molecules investigated using fluorescence microscopy {24 1. A 
similar characterization of the chain shapes at higher semidilute densities appears therefore 
feasible, at least in principle. 
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F. The perimeter length 



As shown in Fig. [TJ 2D chains adopt irregular shapes with perimeters not appearing to be 



smooth, i.e. characterized by a line dimension d p = 1, but clearly fractal (2 > d p > 1) 25 ]. 
In this subsection we analyze quantitatively this visual impression confirming the announced 
key result Eq. ([3]) for the average perimeter length L(s) of chains (s = N) and subchains 
(s < N). 

We define a perimeter monomer as having at least one monomer not belonging to the 
same chain or subchain closer than a reference distance a ~ 1 essentially set by the monomer 
density (see below). Specifically, we have used a = 1.2 in Figs. [TJ |5] and [TD] and for the data 
listed in the Table. The number of such perimeter monomers is called l{s), its mean number 
L(s) = (l(s)) with L(s)/s < 1 being the probability that a monomer of a subchain is on 
its perimeter. (Note that for a continuous chain model a slightly different "box counting" 



method must be used to obtain a finite perimeter length [251].) The main panel of Fig. [TO] 
presents L(s) using the same symbols as in Fig. HI All data collapses on the same master 
curve, confirming nicely the announced exponent 1 — uQ 2 — 5/8 (dashed line) and thus 
a fractal line dimension d p = (5/8) jv = 5/4. This result holds provided that chain and 
subchain lengths are not too small (N, s S> 50). 

Having just confirmed Eq. ([3]) numerically we have still to give a theoretical argument 
to show where it stems from. Using the return probability p r (s) ~ l/s 1+u@2 measured in 
Fig. E] a simple scaling argument can be given following Semenov and Johner js]] . The key 
point is that a monomer in a long subchain cannot "distinguish" if the contact is realized 
through the backfolding of its own subchain or by another subchain of length s. Since the 
probability of such a contact, L(s)/s, must be proportional to p r (s) times the number s of 
monomers in the second subchain we have 

L(s)/s ~ p r {s) x s ~ s- u@2 forl«s<JV (21) 

which is identical to Eq. (]3]). An alternative, but related derivation will be given in Sec- 
tion M£\ 

The fluctuations of the perimeter length are characterized in the inset of Fig. [lTJ We 
present here the histograms P(l, s) of the number of perimeter monomers for different sub- 
chains of lengths s for chains of length iV = 1024. Assuming that the first moment L(s) of 
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the histogram sets the only scale all histograms are successfully brought to a scaling collapse. 
Please note that the fractality of the perimeter does not imply that the histograms have to 
be broad. In fact, they decay rather rapidly as indicated by the two phenomenological 
fits and all moments of the distributions exist. The configurations corresponding to small 
perimeters are rather strongly suppressed (solid line). At variance to this, the decay for 
l(s) — L(s) 3> is found to be successfully described by a Gaussian (dash-dotted line). The 
perimeter fluctuations of different contour sections of these configurations are apparently 
only weakly coupled, if at all. 

The influence of the distance a used to define a perimeter monomer is investigated in 
Fig. [TT] where we present the "relative error" 5L(s)/L(s) = y (l 2 ) / (I) 2 — 1 of the distri- 
butions P(l,s). (The relative error for s = N is listed in the Table.) For clarity, only 
data for one chain length N = 1024 is given for several a as indicated in the figure. Ob- 
viously, if a is too large (a 3> R{s)) all subchain monomers are considered to be perimeter 
monomers, I ~ s, and the perimeter length cannot fluctuate (5L(s)/ L(s) ~ 0). With in- 
creasing s the fluctuations increase first (a > 1.15) and level then off in the limit of large s, 
i.e. SL(s)/L(s) ~ s° in agreement with the scaling found in the inset of Fig. [10J If, on the 
other hand, a is too small (a < 1.15) not all monomers clearly on the contour are detected, 
as can be checked by looking at snapshots similar to Fig. [TJ In this case the fluctuations first 
decay until the subchain length is sufficiently large that the number of detected perimeter 
monomers becomes proportional to the true number. Hence, the relative errors of too small 
and too large a essentially merge for large s or become parallel. The specific value of a is 
thus inessential from the scaling point of view. However, computationally it is important 
to choose a parameter a allowing to probe the asymptotic scaling behavior for as broad an 
s-range as possible. It is for this technical reason that a = 1.2 has been chosen above. 

A method allowing to verify the fractal dimension of the chain contour not requiring such 
an artificial parameter is presented in Fig. [12j We show here the radial pair correlation 
function gi nte r { r , N) between monomers on different chains as a function of the distance r 
between the monomers. The bold line indicates the pair correlation function g(r) between 
all monomers already presented in Fig. [31 The same normalization is used for Winter ( r > N) as 
for g(r), i.e. gi n ter( r , -^0 — > 1 for large distances where both monomers must necessarily stem 
from different chains. Obviously, g(r) > Winter (r, N) for all distances r. For small distances 
r ~ 1 the pair correlation function gi ntcr (r, N) measures the probability that two monomers 
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from different chains are in contact. We remind that for open chains, e.g. self-avoiding 
chains in 3D melts, Winter (f, N) becomes rapidly chain length independent. Our chains are 
compact, however, and only a fraction of the chain monomers is close to its contour. Hence, 
Winter (r, N) must decrease with N for r ~ 1. This is clearly confirmed by the data presented 
in the main panel. From Eq. (I2ip we know already the probability for two monomers from 
different chains to be close to each other, i.e. for both monomers to be close to the chain 
contour. One expects thus to find 

3inter(r, N) ~ L(N)/N ~ AT" 02 ~ iV- 3 / 8 (22) 

for small distances r of the order of a few monomer diameters. This scaling is perfectly 
demonstrated by the data collapse presented for all chain lengths in the inset of Fig. [12] 

G. Intrachain structure factor 

Neither the intrachain size distributions Gi(r, s), nor the bond-bond correlation functions 
-Pi(s) and P2{s) or the contour length L(s) are readily accessible experimentally, at least 
not for classical small-monomer polymers which cannot be visualized directly by means of 
fluorescence microscopy or atomic force microscopy. It is thus important to demonstrate 
that ©2 is measurable in principle from an analysis of the intrachain structure factor F(q) = 
jr J2n m=i( ex P ' ( r « ~~ r m)]) as announced in Eq. @ in the Introduction and as shown 
now in Fig. [13] and Fig. [I4"l 

Figure [12] presents the unsealed structure factors F(q) as a function of the wave vector q 
for a broad range of chain lengths iV as indicated. As one expects, F(q) becomes constant 
for very small wave vectors (F(q) — > N), decreases in an intermediate wave vector regime 
(R g (N) <C 1/q <C 1) and shows finally the non-universal monomer structure for large q 
comparable to the inverse monomeric size ("Bragg peak"). The first striking result of this 
plot is that F(q) does not become chain length independent in the intermediate wave vector 
regime as it does for (uncollapsed) polymer chains in 3D. The second observation to be 
made is that with increasing chain length the decay becomes stronger than the power-law 
exponent —2 indicated by the thin line corresponding to Gaussian chain statistics. 

Since for an open polymer-like aggregate or cluster of inverse fractal dimension v without 
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any sharp surface the structure factor must indeed scale as Q, Q 

F(q) ~ N°q~ 1/u for q > l/R g (N) (23) 

several authors a a i 

111 . Il7j ] have argued that an exponent —2 should be observed for 
2D polymer melts. However, Eq. ( 123]) does not hold for compact structures where strong 
composition fluctuations (of the labeled monomers of the reference chain with respect to 
unlabeled monomers) at a thus well-defined surface or perimeter must dominate the structure 



factor leading to a "generalized Porod scattering" 26) . Since the exponent v = 1/2 for 2D 



melts does not refer to their Gaussian open chain statistics but rather to their compactness, 



Eq. f )23|) is thus inappropriate. Quite generally, the scattering intensity NF(q 
objects is known to be proportional to their "surface" L(N) which implies 26 



28] 



compact 



NF(q) « N 2 / (qR(N)) 2d - dp ~ R(N) d ? ~ L(N) (24) 

where we have used that R(N) ~ N 1 ^ and R(N) dp ~ L(N). For a 2D object with a smooth 
perimeter, i.e. for d p = 1, this corresponds to the well-known Porod scattering F(q) ~ g~ 3 . 
As indicated by the dash-dotted line in Fig. [131 this yields a too strong decay not compatible 
either with our data. Obviously, this is to be expected since we already know from Sec. IIII Fl 
that the perimeter is fractal (d p > 1) and the power-law slope —3 must be a lower bound to 
our data. If we assume, on the contrary, in Eq. ( 1241) a fractal line dimension d p = d — G2, as 
demonstrated analytically in Eq. (121]) . this yields directly the key result Eq. (jl]) anticipated 
in the Introduction. Using = 3/4 we thus have F(q) ~ g -11 / 4 as indicated by the dashed 
line. This power law gives a reasonable fit for the largest chains we have computed. 

The representation of the structure factor used in Fig. [13] is not the best one to check 
the asymptotic power-law exponents and does not allow to verify the iV-scaling implied by 
Eq. ( ]24j) . We have thus replotted our data in Fig. HH using a Kratky representation with 
vertical axis y = (F(q)/N)Q 2 and a reduced wave vector Q = qR g (N). Using the measured 
radius of gyration R g (N) given in the Table this obviously allows to collapse all data in the 



Guinier regime for Q ^ 1 where 36] 



2/(Q) = Q 2 (l-?). (2r,) 



d 

The observed data collapse is, however, much broader in Q and the more the larger the chain 
length. The deviations observed for large Q are due to (chain length independent) physics 
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on scales corresponding to the monomer size ( "Bragg peak" ) already seen in Fig. [T3j The 



Debye formula for Gaussian chains |36j is given by the thin line which becomes constant in 
this representation for Q 3> 1 in agreement with Eq. ( I23p . At variance to this, a striking 
decay of y(Q) is observed over a decade in Q confirming observations by Yethiraj [11] and 
Cavallo et al. [15;] using much shorter chains. Due to the scaling of y(Q) as a function 
of Q this decay implies the chain length dependence seen for the unsealed structure factor 
in Fig. [13j With increasing chain length our data approaches systematically the power 
law y(Q) ~ Q~® 2 given by Eq. (jlj) and indicated by the dashed line. Even longer chains 
obviously are warranted to unambiguously show the predicted asymptotic exponent —3/4 
in a computer experiment. 

In the preceding two paragraphs we have used the fractal line dimension d p = d — 2 
derived via Eq. fl2~Tj) together with the generalized Porod scattering scaling Eq. f l24|) to 
demonstrate the key result Eq. 01]). Interestingly, the structure factor can be computed 
directly without the scaling argument Eq. ( l2T|) using that G e (r, s) ~ G 2 (r, s) for s <C N as 
discussed at the end of Sec. IIIIBI For asymptotically long chains the structure factor thus 
can be well approximated as 

i r N 

J ds2(N-s)G 2 (q,s) (26) 

using the Fourier transform G 2 (q, s) of G 2 (r, s). Within the Redner-des Cloizeaux approxi- 
mation, Eq. OH]), this yields an analytic formula, Eq. QA12I) . given in the Appendix. Readily 
computed numerically, this theoretical prediction is represented by the solid line in Fig. [TU 
Since Eq. (JHJ) and Eq. Q26l) are both approximations this result is not strictly rigorous. 
However, by construction our formula must yield the Guinier regime, Eq. fl25|) . for small 
Q = qR g (N) and since for large Q only the 02-exponent matters for large iV it is only 
around the hump Q w 2 where deviations could be relevant. Fig. [141 shows that in practice 
our approximation agrees well for all Q as long as the wave vector q does not probe local 
physics (Bragg regime). 

As shown in the Appendix, the Redner-des Cloizeaux approximation Eq. ( 1A12I) reduces 
to a power law for wave vectors Q ^> 1 corresponding to the power-law regime of Eq. ([7]), 

2 / 3 \ -( 1+Q 2/ 2 ) i qs 

*«> K w^m (jtW q " W- (27) 

in agreement with the key Eq. (j3J) given in the Introduction. Equation (I27p is represented by 
the dashed lines in Fig. (jl3~|) and Fig. QT4|) . Comparing Eq. Q24"|) with Eq. Q27|) demonstrates 
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that 2D melts are characterized by a fractal line dimension 

d p = d - G 2 = 5/4. (28) 

Hence, using a slightly more physical route as the scaling argument given in Sec. IIIIFI we 
have confirmed the fractal line dimension of the chain perimeter L(N). By labeling only 
the monomers of sub-chains (which corresponds to a scattering amplitude sF(q) ~ L(s) ~ 
R(s) dp ) the above argument is readily generalized to the perimeter length L(s) of arbitrary 
segments of length s < N. 

IV. CONCLUSION 

Using scaling arguments and molecular dynamics simulation of a well-known model 
Hamiltonian we investigated various static properties of linear polymer melts in two di- 
mensions. We have shown that the chains adopt compact conformations [y = 1/d = 1/2). 
Due to the segregation of the chains the ratio of end-to-end distance R e {N) and gyration 
radius R g (N) becomes R%(N)/R^(N) ~ 5.3 < 6 (Fig. ED and the chains are more spherical 
than Gaussian phantom chains (Fig. [HD- More importantly, it is shown that the irregu- 
lar chain contours can be characterized by a fractal line dimension d p = d — Q 2 = 5/4 
(Figs. [TUl and [T2|) . This key result has been demonstrated analytically using two different 
scaling arguments given in Sec. IIIIFI and Sec. IIII Gl both based on the numerically tested 
power-law scaling of the intrachain size distribution G c (r,s) Gz(r, s) ~ r 02 for small 
distances r <C R e (s) with G 2 = 3/4 (Fig. ED- Compactness and perimeter fractality repeat 
for subchains of arc-lengths s down to a few monomers due to the self-similar structure of 
the chains (Figs. |2|41 and [TO]) . Measuring directly the return probability of the chain after s 
steps, the second Legendre polynomial ^(s) of the bond vectors decays as -P 2 (s) ~ l/s 1+u@ ' 2 
(Fig. [HD- Interestingly, as implied by the generalized Porod scattering of a compact object 
with fractal "surface", Eq. (l24j) . the predicted fractal line dimension should in principle be 
accessible experimentally from the power-law scaling, Eq. (j3J), of the intrachain structure 
factor F(q) in the intermediate wave vector regime. Computationally very demanding sys- 
tems with chain lengths up to iV = 2048 have been required to test the proposed scaling of 
the structure factor (Fig. [T41) . 

We would like to stress that our results are not restricted to a particular melt density, 
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but should also hold for all densities provided that the chains are sufficiently long to allow a 
renormalization of all length scales in terms of semidilute blobs This is of some interest 
since chain conformations of semidilute 2D solutions of large-monomer polymers (such as 



DNA or brushlike polymers) are experimentally better accessible than dense melts [17|, |20 |. 
Obviously, these macromolecules are rather rigid and in view of the typical molar masses 
currently used, deviations are to be expected from the asymptotic chain length behavior we 



11 



12 



191 ] it should thus be 



focused on. Following previous computational work jf], |8. 9. 
rewarding to reinvestigate the scaling of flexible and semiflexible chains in 2D semidilute 
solutions to see how finite- N effects may systematically be taken into account. 

Interestingly, the fractality of the perimeter precludes a finite line tension and the shape 
fluctuations of the segments are not suppressed exponentially [6], but may occur by ad- 
vancing and retracting "lobes" in an "amoeba-like" fashion. This opens the possibility for a 
relaxation mechanism, specific to 2D polymer melts, in which energy is dissipated by friction 
at the boundary between subchains. Following a suggestion made recently |3|, the longest 
relaxation time r(s) of a chain segment should, hence, scale as 



t(s) ~ L(s) 3 ~ s a with a = 3(1 - uQ 2 ) = 15/* 



(29) 



rather than with a = 2 as predicted by the Rouse model which is based on Gaussian chain 
statistics 36(. As Gaussian chain statistics is inappropriate to describe conformational 
properties of 2D melts, there is no reason why a modeling approach based on this statistics 
may allow to describe, e.g., the composition fluctuations at the chain contour. Since the 
latter can in principle be probed experimentally using the dynamical intrachain structure 
factor F(q,t) 26( this is an important issue we are currently investigating. 
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Appendix A: Calculation of intrachain structure factor 

The intramolecular structure factor F(q) may be rewritten generally as [3j 

F(q) = i J ds 2(N - s)G e (q, s) (Al) 

using the Fourier transform G e (q,s) of the two-point intramolecular correlation function 
G e (r,s) averaging over all pairs of monomers (n,m = n + s — 1) discussed in Sec. IIIIB1 
As we have seen in Sec. IIII B[ G e (r, s) is well approximated by the distribution G2{r, s) for 
s <C N. For asymptotically long chains it is justified to neglect chain-end effects (s — > N), 
i.e. physics described by the contact exponents O and ©i. Assuming thus translational 
invariance along the chain contour the structure factor is given approximately by 

i r N 

F(q) ~ ^ j ds 2(N - s)G 2 (q, s), (A2) 

the factor 2(N — s) counting the number of equivalent monomer pairs separated by an arc- 
length s. Using the Redner-des Cloizeaux approximation, Eq. (jSJ, for i = 2 we compute 
first the 2D Fourier transform 

POO 

G 2 (q,s)= c 2 x @2 e- k2x2 2irxdxJ (qx) (A3) 

I — I 
with 2nJ (z) = 2 J °° cos(z cos(9))d9 being an integer Bessel function |4l| and x = r/R 2 (s) = 

r/b c s 1 ' 2 , 6 2 = 3/4, k 2 = 1 + 6 2 /2, c 2 = k k 2 2 /nT(k 2 ) as already defined in Sec. ImBl As 



can be seen from Eq. (11.4.28) of Ref. 4l|, this integral is given by a standard confluent 



hypergeometric function, the Kummer function M(a,b, —z), 

G 2 (q,s) = M(l + e 2 /2,l,-z) (A4) 



with z = q 2 b1s/Ak 2 . According to Eq. (13.1.2) and Eq. (13.1.5) of [4jJ the Kummer function 
can be expanded as 

OjZ 

M(a, b,-z) w 1 - — for \z\ < 1, (A5) 
b 

M{a, b, -z) w -^—z- a for z > 1. (A6) 
T(b — a) 

Using Eq. (1A4I) this yields, respectively, the small and the large wave vector asymptotic 
behavior of the Fourier transform of G 2 (r, s) 

G 2 (q, s) w 1 - (1 + 2 /2)^ for z < 1, (A7) 
G 2 (g, 5) « T ^^r^ (1+02/2) ~ <T (2+02) for « » 1. (A8) 
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Note that Eq. (1A7j) implies G 2 {q = 0, s) = 1 as one expects due to the normalization of 
G 2 (r,s). 

After integrating over s following Eq. flA2|) and defining Z = q 2 b 2 N/Ak 2 one obtains for 
the Guinier regime of the structure factor 

F(q) « N ( 1 - 1 + f 2 ^ 2 Z) for Z < 1, (A9) 



3 

i.e. according to Eq. ( )25l) we have, as one expects, 

Rl(N) = >i±f^ = (MO, 

This is consistent with Eq. f[T2"j) and the ratio (b e /b g ) 2 = 6 with 6 e determined from subchains 
with s <C N, Eq. pip . Eq. flAlOj) is of course slightly at variance with the measured ratio 



(Re(N) I R g (N)) 2 = (b e +/b g ) 2 < 6 due the end-effects not taken into account in Eq. (1A2I) . 
The power law behavior of the structure factor for large wave vectors announced in Eq. (j3J) 
is obtained by integrating Eq. (IA80 with respect to s. This gives 



F (q) « - 2N z-^l 2 ) ~ N -^/2 q -(2+e 2 ) for z ^ L (An) 

r(2 — 02/2) 

Obviously, it is also possible to directly integrate Eq. ( 1A4j) with respect to s according to 
Eq. (1A2I) . This yields the complete Redner-des Cloizeaux approximation of the structure 
factor 

F <«' - 2M(l + %,2,-Z 



\ 2 
- m(i + ^,3,-Z 

+ 1( 1 + y) z m ( 2 + T' 4 ' (A12) 

which can be computed numerically. Using again the expansions of the Kummer function, 
Eq. (1A5h and Eq. (1A6I) . one verifies readily that Eq. ( 1A12I) yields the asymptotics for small 
and large wave vectors already given above. It is convenient from the scaling point of view 
to replace the variable Z used above by the reduced wave vector Q = qR g (N) substituting 

Z^Q 2 (l + 6 2 /2) = -Q 2 , (A13) 



as suggested by Eq. (lAlOj) . This gives the curve represented by the solid line in Fig. [T 



Eq. (lAllj) reexpressed in these terms is given by Eq. (T2"T|) in the main text. Note that 
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due to this substitution the Guinier limit of the Redner-des Cloizeaux approximation of the 
structure factor is correct by construction. 
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N 


R c 


Rg 


(f) 2 


<Ai> 

\ A 2/ 


Ai 


A 2 


L(N) 

N 


SL(N) 
M iv ) 


32 


8.1 


3.4 


5.7 


4.9 


0.66 


0.56 


0.55 


0.21 


64 


11.7 


5.0 


5.4 


4.7 


0.65 


0.54 


0.44 


0.22 


128 


16.7 


7.2 


5.4 


4.6 


0.64 


0.54 


0.35 


0.23 


256 


23.8 


10.3 


5.3 


4.5 


0.64 


0.53 


0.27 


0.23 


512 


34.0 


14.7 


5.3 


4.5 


0.64 


0.53 


0.21 


0.23 


1024 


48.2 


20.8 


5.3 


4.5 


0.63 


0.52 


0.16 


0.23 


2048 


66.4 


28.9 


5.3 


4.4 


0.62 


0.51 


0.13 


0.23 



TABLE I: Various conformational properties defined in the main text for 2D polymers at monomer 
number density p = 7/8 as a function of chain length N. The aspherity of the chains is characterized 
by the aspect ratio (Ai) / (A 2 ) and the moments Aj and A 2 of the eigenvalues Ai and A 2 of the 
inertia tensor. Note that Rg = (Ai) + (A 2 ). The mean perimeter length L(N) and its relative 
fluctuation 5L(N)/L(N) given in the two last columns have been obtained assuming a reference 
distance a = 1.2 (as indicated in Fig. [3]) within which at least one monomer from another chain is 
to be found to identify a monomer as a perimeter monomer. 
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FIG. 1: Snapshot of a 2D polymer melt with chain length N = 2048 and monomer number density 
p = 7/8. Due to the excluded volume interactions the monomers do barely overlap and chain 
intersections are strictly forbidden. One chain in the middle is fully drawn while for the other 
chains only the perimeter monomers interacting with other chains are indicated. A perimeter 
monomer is defined here as having at least one monomer not belonging to the same chain within a 
distance r < a = 1.2. Numbers refer to an arbitrary chain index used for computational purposes. 
The chains are compact, i.e., they fill space densely, and interact typically with about 6 other 
chains. However, compactness does apparently not imply a disk-like shape which would minimize 
the perimeter length L(N). Quantitative analysis reveals that the contours are fractal with a 
fractal line dimension d p = 5/4. 
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N=2048,s=256 



chain "1" 




chain perimeter 



FIG. 2: Self-similarity of compactness and perimeter fractality on all scales shown for the chain 
"1" from the snapshot Fig. [TJ The solid line indicates the perimeter of this chain with respect to 
monomers of other chains. We consider 8 consecutive subchains of length s = 256 and compute 
their respective perimeter monomers being close to monomers from other chains or subchains. The 
subchains are compact and of irregular shape, just as the total chains (s = N) presented in Fig. [TJ 
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FIG. 3: Radial pair correlation function g(r) between all monomers for chains of length N = 2048 
at monomer number density p = 7/8. Note that g(r) — > 1 for large monomer distances r. The 
correlation function oscillates strongly with a period given by the bead diameter a = 1. As shown 
in the inset, the amplitude of these oscillations decays exponentially as expected for a dense liquid 
without long-range positional order [34]. The arrow indicates the position of the reference length 
a = 1.2 used for identifying a monomer as being in contact with a monomer from another chain or 
sub chain. 
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FIG. 4: End-to-end distance R e (s) = 
of "subchains" or "chain segments" 



(r 



2\l/2 



(top data) and radius of gyration R g (s) (bottom data) 



371 ] containing s = m — n + 1 monomers as indicated by the 
sketch on the left hand side. The data is averaged over all s-subchains possible in a chain of length 
N . Open symbols refer to subchains with s < N for N = 1024 (squares) and N = 2048 (spheres), 
stars to overall chain properties (s = N). The indicated power-law slopes confirm Eq. (j2J), i.e. the 
chains are compact on all scales. 
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(a) G (r,N) G,(r,W2) G 2 (r,N/2) G e (r,s) 





=n+s-1 



9 = 3/8 1= 1/2 9 2 = 3/4 



s « N 



e 2 = 3/4 




FIG. 5: (a) Sketch of the intrachain probability distributions Gj(r, s) described in the main text and 
the contact exponents 0« predicted by Duplantier [2]. (b) Scaling plots are presented for Go(r, N), 
Gi(r,N/2) and G 2 {r,N/2) for N = 1024 (squares) and N = 2048 (spheres) and for G c (r,s) for 
N = 1024 3> s with s = 256 (triangles) and s = 512 (diamonds). All data for different N and 
s collapse on the respective universal master curves if y = RfGi(r,s) is plotted vs. the reduced 
distance x = r/Rj, with i2? being the second moment of the corresponding distribution. The thin 
line indicates the Gaussian distribution y = exp(— x 2 )/ir expected for ideal chains in 2D. The power 
laws y ~ x 0i observed for x <C 1 (dashed lines) confirm Duplantier's prediction. The solid line at 
the bottom shows the Redner-des Cloizeaux formula, Eq. (|8|), for G2(r,s = N/2) ~ G e (r,s). 
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FIG. 6: Reduced typical subchain size R(s)/s 1 ^ 2 as characterized by the root-mean-square end- 
to-end distance R e (s) (top data) and the radius of gyration R s (s) (bottom data) using the same 
symbols as in Fig. 21 Interestingly, R e (s)/s 1 ^ 2 decays for s > N/2 due to chain-end effects. The 
solid horizontal lines indicate the effective segment sizes b e = 1.6 and b g = 0.65 obtained for 
100 <C s <C N from R e {s)/s 1 / 2 and R g (s)/s 1 ^ 2 , respectively. The dashed line corresponds to an 
apparent effective segment size 6 e * = 1-5 from the end-to-end distances R e (N) of our longest chains. 
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FIG. 7: Bond-bond correlation function P\(s) = (e n ■ e m ) vs. curvilinear distance s — 1. Main 
panel: The first Legendre polynomial shows an anti-correlation at s — 1 ~ 13 due to the local 
wrapping of the chain. Inset: A second anti-correlation peak is visible at s ~ N/2 which is caused 
by the reflection of the confined chain. Due to Eq. (|14p this peak is related to the decay of the 
reduced subchain size R^{s)/s 1 / 2 visible in Fig. [6j 
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FIG. 8: Second Legendre polynomial P2 = ((e n • e m ) ) — 1/2 as a function of curvilinear distance 
s — 1 (a) and as a function of intrachain distance r (b). (a) P2(s) decays over two orders in 
magnitude as a power law (dashed line) with an exponent 1 + U62 = 11/8 in agreement with 
a return probability given by Eq. f) 1 5 1) . (b) 2irrP2{r) oscillates strongly with an exponentially 
decreasing amplitude. Being dominated by a nematic peak at the origin (r « 1) it acts as a 
5-function. 
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FIG. 9: Aspherity of chains (stars) and subchains for N = 1024 (squares) and N = 2048 (spheres) 
using log-linear coordinates. The aspherity is characterized by the moments Ai and A2 of the 
eigenvalues Ai and A2 of the inertia tensor discussed in the main text. Both moments decay weakly 
with s and N. The horizontal lines are guides to the eye indicating possible plateau values for 
asymptotically long chains. 
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io 1 io 2 s io 3 

FIG. 10: Mean perimeter length L(s) = {I) (main panel) and perimeter length distribution P(l, s) 
(inset) for a = 1.2. In agreement with Eq. (J3]) the perimeter scales as L(s) ~ s 1 ~ u ®2 ^ s 5/8 ( do ^ 
line). Stars refer to perimeters of chains (s = N), open symbols to subchains of lengths s < N 
in chains of length N = 1024 (squares) and N = 2048 (spheres). The histograms for different s 
collapse if rescaled using L(s) as presented in the inset for different segment sizes and iV = 1024. 
The distributions are lopsided decaying much sharper for small-perimeter segments (l/L(s) -C 1). 
The exponential cut-offs indicated for both limits are phenomenological fits. 
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FIG. 11: Relative error 8L{s) / L{s) as a function of s for chains of length N = 1024 and for different 
distances a as indicated. The relative error increases for small s and too large a-values. It decreases 
for a < 1.15 since these values are too small for the given density and not all perimeter monomers 
are detected. The relative error becomes constant for large s irrespective of the value of a. Due to 
chain end effects the relative error increases slightly for s ~ N. 
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FIG. 12: Radial pair correlation function Winter ( r , N) between monomers on different chains using 
the same symbols as in Fig. [8] for the chain length N. Also indicated is the pair correlation function 
g(r) between all monomers for N = 2048 (solid line). Measuring the probability that a monomer is 
close to the chain perimeter Winter ( r ) N) decreases strongly with N. This decay is perfectly described 
by <?mter( r , N) ~ _/v~ 1/6 ' 2 = as shown by the scaling collapse of the data presented in the 

inset. 
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FIG. 13: Intramolecular structure factor F(q) as a function of wave vector q for different chain 
lengths ./V as indicated. In striking contrast to 3D melts F(q) becomes chain length independent 
only for wave vectors corresponding to the monomer scale ("Bragg peak"). To characterize the 
decay in the intermediate wave vector regime our data is compared with three power-law exponents 
—2, —11/4 and —3 indicated by the thin, the bold dashed line and the dash-dotted line, respectively. 
The first exponent corresponds to Eq. (j23|) . the second exponent to Eq. (j3J) and the last exponent to 
Eq. (j24|) . i.e. to the Porod scattering of a compact 2D object with smooth surface (d = 2,d p = 1). 
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Q = qR g (N) 

FIG. 14: Kratky representation of the structure factor F(q) tracing y = (F(q)/N)Q 2 as a function 
of the reduced wave vector Q = qR g (N) for different N using the same symbols as in Fig. [131 The 
Debye formula (top thin line) corresponds to a plateau for Q > 1. At variance to this a strong 
non-monotonous behavior is revealed by our data which approaches with increasing N a power law 
exponent —02 = —3/4 (dashed line) corresponding to a compact object of fractal line dimension 
d p = d — 82 = 5/4. The Porod scattering for a compact 2D object with smooth perimeter is given 
by the dash-dotted line, the Fourier transform of the Redner-des Cloizeaux approximation by the 
solid line. 
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